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| These lectures discuss aspects of high energy evolution in QCD. This 

£T} , includes the derivation of the JIMWLK equation, basic physics of its solu- 

£N| ' tions and recent work on inclusion of Pomeron loops. The entire discussion 

is given in the Hamiltonian framework which gives direct access to the 

■ evolution of hadronic wave function under Lorentz boost. 
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', 1. Introduction 

These lectures deal with the description of hadronic scattering at high 
Oh! energies, and more specifically with the evolution of the scattering matrix 

with energy pP,[5],jS]- It has been a very lively and active field for the last 10 
years [3], 0)0 and especially lately there has been a new burst of activity. 
The latest developments are related to what is sometimes called Pomeron 
^ ■ loopsjHj, or Fluctuations 9 or alternatively Saturation effects in the projec- 

tile wave function jSJ. It is the last characterization that is most suited to 
the nature of these lectures. Somehow, although the physical underpinning 
of the low x (or high energy) evolution is the understanding how hadronic 
wave functions evolve with energy, the language of the wave function is not 
frequently used in this context. I believe that one can learn a lot by dealing 
directly with the wave function, since it contains all the nontrivial informa- 
tion there is to have, so my heart goes out to the wave function. I have 
therefore decided for these lectures to adhere to its point of view on Life, 
the Universe and Everything, and as a small part of Everything, to show 
how to think about the evolution and derive the evolution equation from 
this perspective. 

The high energy evolution in its latest incarnation became more or less 
synonymous with saturation. Why is it interesting to study high energy 
evolution? One obvious answer is of course that there is an experimental 
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drive to understand the high energy regime. The resurgence of interest in 
this topic in the 90's was firmly triggered by the DIS data from HERA, which 
extended experimental reach to x ~ several x 10 -5 . The striking observation 
there was that the at low x (high energy) the number of gluons in the proton 
shoots up very fast (about 20 gluons at Q 2 = 20 GeV, x = 10 -4 ). How and 
why does this number grow so fast? There is still no definitive answer to 
this question, in the sense that there are more than one fits to HERA data, 
some based on saturation^! and some not^J, and it is strictly speaking not 
possible to differentiate between the physical mechanisms involved on this 
basis. Nowadays this motivation is still there and is perhaps even stronger, 
since LHC is hopefully coming on line soon. There much larger numbers of 
gluons will be created in energetic p-p collisions. Another motivation comes 
from RHIC. Although the energies there are not extremely large, due to 
the fact that the colliding objects are nuclei the effects of the evolution are 
expected to appear much sooner. It may well be that large gluon densities 
are at work there. Some features of the data are at least qualitatively 
described by the saturation models \TQ, which are based in some aspects on 
low x evolution|15j. And finally, Pb-Pb collisions at LHC will lead to even 
higher gluon densities. 

True as it may be, at the moment the high energy evolution and satura- 
tion physics can not claim firmly that it is of phenomenological relevance. 
For me personally the more appealing side of this field is that the physics 
is of fundamental theoretical interest. 

We all know that QCD is nonperturbative in the Infrared. Masses, 
decay widths, elastic scattering cross sections can not be obtained without 
recourse to nonperturbative physics. But in the Ultraviolet the theory is 
perturbative. Quarks and gluons when probed on a short distance and time 
scales behave as if they are almost free and thus some physical processes can 
be calculated perturbatively. The prime example of this is DIS. Electron 
scatters off a nucleon with large (space like) momentum transfer q 2 = —Q 2 . 
A simple picture of this process is given by the parton model. The proton 
comes in with momentum in the z direction and is hit by a space like 
photon (emitted off the incoming electron). By measuring the outgoing 
electron one has access to two kinematical variables, the energy transfer 
from the electron to the hadronic system, and the scattering angle of the 
electron. The two Lorentz scalar quantities which encode this information 

are Q 2 and x = ^p.q ■ The parton model is formulated in the infinite 
momentum frame of the proton, where the four momenta of the proton and 
the virtual photon involved in the collision are 



P IM ( P ) = (E p = P,0,0,P), P^cx) 



zakopane printed on February 2, 2008 



3 



e 



Y 



P 




X 



Fig. 1. Deeply Inelastic Scattering 

M7) = (^ = |^,Q[i + J^] 1/2 ,o) 

The energy of the colliding 7 — p system is s = (Q^ + P^) 2 = Q 2 (- — 1) - at 

small x obviously s = hence small x limit is the same as large energy 
limit in this case. 

The physical picture of the parton model is simple. In the frame we 
have chosen, the photon is a very small probe of transverse size ~ ^. What 
does it see when it hits the proton? There is of course a lot of stuff inside 
the proton, which is mainly soft muck. But a small probe does not see this 
- it whizzes through the muck unless it strikes something hard. That is 
to say, it only scatters off the proton if it encounters another small object 
of roughly the same size as itself. The object in this case of course must 
be a quark, since only quarks are charged and interact with the photon. 
The scattering event itself therefore is a short distance process which also 
happens very quickly in time, with the relevant distance-time scale of order 
l/Q. At large Q 2 this is perturbative. 

This short distance processes are described perturbatively, but such hard 
fluctuations over the soft background are very rare, thus the cross section 
is small. In other words there is only a small number of quarks in the 
proton that a hard probe can see. The cross section for the 7 — p subprocess 
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(forgetting about the 7 - emission vertex from the scattered electron) is easy 
to write down 



7 scatters if it hits the quark, and N g is the average number of the hard 
quarks in the proton. The parton model interpretation of the variable x 
is the longitudinal momentum of the struck quark, or more precisely the 
fraction of total proton momentum that the struck quark happens to carry 
- hence its place of honor as a variable of which N depends. The source of 
this interpretation is the kinematics of the two body 7 — 5 collision. Think 
of the struck quark to be on shell, massless particle. Before the collision 
it carries momentum (k, 0, k) since it moves in the same direction as the 
proton. After the collision its four momentum is k! 



But (k'f = 0, so (k + ^pf - Q 2 - k 2 = - Q 2 = 0. Thus k = xP. 



The unknown quantity in the cross section is of course the number of 
quarks. We do not know it, since it is a nonperturbative quantity. It is some 
property of the proton wave function, which develops over a long period of 
time and thus crucially depends on the nonperturbative physics. However 
one can certainly ask within the perturbation theory how it changes, when 
we change the "resolution scale" Q, as long as Q is large. If Q increases, the 
change in the quark number is perturbative - it changes simply due to very 
fast fluctuations. In a way, the old quarks and gluons "split" by perturbative 
radiation processes. As we increase the resolution of the probe, and also 
decrease the time duration of the interaction, the probe resolves more and 
more of these fast fluctuations in the wave function. This is described by 
the QCD DGLAP evolution equation. In this evolution the picture of the 
hadron basically does not change. There are a little bit more partons as we 
go to higher resolution (the growth is logarithmic in perturbative QCD), 
since when we look closer sometimes we see that what we thought was one 
parton, is in fact two sitting close together. But the density of partons in 
the transverse space actually decreases, since the size of the resolved partons 
decreases much faster than their number increases. 

But what happens if we increase the energy of the process keeping the 
resolution fixed (decrease x at fixed Q 2 )l Again there will be a growth of 
number. The partons split and multiply by the same process of radiation of 
gluons. But the splitting now is predominantly in the longitudinal direction 
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Fig. 2. Cartoon of the DGLAP evolution. With improved transverse resolution the 
number of partons grows but the density decreases. 

and on the same transverse scale. The partons split, but keep their trans- 
verse size. When looked "face on" the proton becomes denser and denser. 
Again this is the effect of fluctuations. This time as it turns out (and as we 




Fig. 3. Cartoon of the low x evolution. As the energy increases the partons of hxed 
transverse size multiply. The density in the transverse plain grows. 

will see later) the dominant contribution comes from gluons. So think about 
it as if "valence" partons are surrounded by gluon cloud. The gluons in this 
cloud fluctuate very quickly, and are not seen by the probe, since their effect 
averages to zero during the time of the interaction. But increasing energy 
means that we are boosting our proton. Boost leads to Lorentz time dila- 
tion, so what used to be fast becomes slow. The gluons "freeze" . This Cold 
Gluon Cloud (CGC) 1 is indeed seen by the probe and is the main culprit for 
the growth of cross sections with energy. The caveat is of course, that the 
photon does not see gluons directly, and so for DIS this effect comes into 

1 You have probably noticed that the abbreviation for Cold Gluon Cloud - CGC is 
the same as for Color Glass Condensate introduced in This is in fact fortunate, 
since both refer to exactly the same thing. I personally prefer Cold Gluon Cloud as 
it expresses better the essence of the relevant physics. 
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being slowly - only due to the fact that the sea quark distributions grow due 
to the growth of gluon cloud. But for p — p collisions for example one has 
direct gluon-gluon subprocesses, and those trigger directly on the gluons in 
the proton. 

Evolving far enough in energy will create a system where many gluons 
overlap in the transverse plane. The physics of the proton wave function 
becomes quite different - physics of a dense gluonic system, rather than that 
of a dilute small number of partons. One has to be a little careful with the 
terminology - we are still looking at this system on relatively short time 
scales, on which the gluon cloud is static. Thus even though the system is 
dense, the physics is not that of thermodynamic equilibrium. Nevertheless 
the high density certainly has profound effect on all kinds of observables. 

What can we expect to see, if we collide two such dense "black" objects? 
Think of the limit when the whole face of the proton is covered by these 
overlapping frozen gluons. It is more or less homogeneous on the low mo- 
mentum scale, but is still of course grainy if we probe it at short enough 
distance. One can define the scale, which is proportional to the transverse 
density of gluons. Let's call it the saturation momentum - Q s . This scale 
separates the "homogeneous" regime from the "grainy" regime. For spatial 
resolution worse than 1/Q S gluons strongly overlap in the transverse plane, 
while at resolution better than 1/Q S the gluonic system looks dilute. So 
l/Qs is the resolution at which the packing fraction of gluons in the trans- 
verse plane is about one half. As this is the characteristic transverse scale 
in the problem, we expect that after collision most final state gluons will 
pick up this momentum k± ~ Q s . If Q s is large the physics at this scale 
is perturbative - we produce hard "minijets". These minijets will dominate 
the final state and so maybe we will even be able to describe the bulk of the 
cross section (total cross section?) in terms of these perturbative processes. 
They are certainly not trivially calculable, as we have to learn how to deal 
with dense gluonic medium. Clearly naive perturbation theory in the sense 
I described above does not work in this dense medium. But at least in prin- 
ciple we do not require the knowledge of soft physics of confinement and we 
expect with some effort to be able to figure the problem out. 

Of course, the k± ~ Q s is not the only type of physics that exists 
even in this regime. There are even more perturbative processes, like real 
hard jets. Those will stay perturbative in the usual sense, although their 
multiplicities will be also affected, as presumably gluon distributions will 
change significantly in this regime. There will also remain the soft physics. 
Proton has finite size, and clearly the density will grow in the center faster 
than on periphery. There will always be a peripheral region which is "gray" 
and not "black" , that is where the analog of Q s is of order Aqc d ■ 

We cannot hope to be able to describe the peripheral gray region in any 
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Fig. 4. Schematic portrait of a hadron at high energy. A black central region is 
surrounded by "gray" disc where the saturation momentum Q s is of order Aqcd- 

kind of perturbative approach. Confinement rules there, and confinement 
therefore will rule eventually for total and elastic cross sections. However 
it is conceivable that in some intermediate energy regime, the bulk of the 
proton is already "black" , but the peripheral region is still not too large. In 
this regime the total cross section can be dominated by the minijet processes, 
and might be describable perturbatively! It would be almost too good to be 
true if this did happen. This is a very long shot, and we would be extremely 
lucky. But in principle such a possibility exists and is very interesting to 
explore. 

To summarize, some of the questions we may hope to address in the 
high energy - large density regime are 

1. How does this blackness develop in the first place (Pomeron loops are 
presumably important)? 

2. How does the density of the gluon cloud change with energy, and 
what is the density profile (transverse momentum dependence)? 

3. What is the final state structure when colliding these "black" objects. 

4. What can we say about the total cross section? 

I will certainly not answer all of these questions in these lectures. In 
fact I will not even address most, and for those I will address only tentative 
answers will be provided. Nevertheless a lot has been done and some has 
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been understood in this area and there is ample room for further progress, 
which is enough of a motivation in itself. 

So let's start climbing into the Cold Gluon Cloud. We all fly a lot these 
days and know that going through a cloud could be often bumpy. So fasten 
your seatbelts and let's go. 

The plan of these lectures is the following. 

a) Derive the evolution equation (evolution in energy, x or rapidity) for 
the scattering of a dilute perturbative projectile on a dense target. This 
is what is variably known as BK, JIMWLK, B-JIMWLK and B-JIMWLK 2 
equation. 

b) Discuss briefly what is known about the qualitative behavior of its 
solutions. 

c) . Discuss what kind of corrections have to be taken into account 
when the projectile itself becomes dense (Pomeron loops). Calculate some 
corrections of this type and discuss a little bit what kind of physics we 
expect from them. 

d) . Prove an interesting duality relation which the complete evolution 
equation, including the Pomeron loops must satisfy. 

2. JIMWLK as we know it. 

Consider a small perturbative projectile with the wave function \P > 
impinging from the left on a large dense target with the wave function \T >. 
The target is dense, and therefore has large gauge fields. The projectile on 
the other hand is just a bunch of partons. In the following I will disregard 
quarks, and will only consider gluons as the dynamical degrees of freedom. 

We assume that both the target and the projectile move with close to 
velocity of light. Their longitudinal sizes are therefore Lorentz contracted 
in the lab frame into pancakes of almost zero width. Thus all the projectile 
partons are localized around x~ = 0, while the color fields of the target are 
at x + = 0. The only nontrivial dynamics therefore happens in the transverse 
plane. 

We want first of all to derive an expression for the S'-matrix for the 
collision of two such objects, and then see how it evolves with energy. We 
assume that the initial energy is already quite large. 

2.1. The forward amplitude 

First consider a single gluon in the projectile wave function. It has a 
large energy and therefore while traveling through the target it continues 
along the straight line. Of course it could happen from time to time that it 
experiences hard scatterings and changes its direction of motion, but these 
events are rare and contribute very little to the total cross section. So the 
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single most important thing that happens to the traveling gluon as a result 
of the interaction is that its color degree of freedom rotates. 

\x,a)=af(x)\0)^S ab (x)\x,b) (2) 

The 5-matrix operator in this approximation is block diagonal. It is di- 
agonal in the space of transverse coordinates, particle number etc, and is 
nontrivial only in the space of color indices. The matrix S ab must be unitary, 
since the ^-matrix is a unitary operator. The forward scattering amplitude 
is defined as the overlap of the \in) and \out) states, and thus is given by 
the unitary matrix we have just defined 

in (gluon, 6|gluon, a) out = (0|<(x)5af(x)|0) = S ab (x) (3) 

Now in fact one can express the matrix S in terms of the target fields 
through which the gluon is propagating. The target of course is a quantum 
object, and therefore the fields in it fluctuate. However the propagation 
happens very fast, as we assume that the energy of the projectile is large. 
Therefore during the propagation the projectile really sees a fixed field con- 
figuration. Since we are considering the projectile to consist of partons, it 
is convenient to use the light cone gauge, in which the wave function of 
the projectile has partonic description, namely A + = 0. In this gauge nat- 
urally the largest component of the field in the target is A~ . The target 
moves with large rapidity in the negative 23 direction. Recall that under 
the Lorentz boost in the negative 23 direction the components of vector 
potential transform as 

A~ [x~ , x + , x±) — ► e Y A~ (e Y x + , e~ Y x~, x±) 

A + (x~, x + , x±) — ► e~ Y A + (e Y x + , e~ Y x~ , x j_) 

Ai{x~ , x + , x±) — > Ai(e Y x + , e~ Y x~ , x±) (4) 



where Y is the boost parameter. Thus a large A~ field component is gen- 
erated in the target when accelerating it. The S matrix of a fast particle 
interacting with the A~ field is given by the eikonal formula 



S ab (x) = 



PexpS^ig J dx + a c (x,x + )T c 



ah 



(5) 



where g is the strong interaction coupling constant, a = A~ (change of 
notation for historical reasons only), T£ c = if abc is the SU(N) generator 
in the adjoint representation and P denotes the path ordering along the 
trajectory of the propagating gluon.. Now, the fields a(x) are of course 
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operators in the target Hilbert space. So to calculate the 5-matrix for the 
process, we have to average over the target wave function. 

a(l gluon) = (T\S(a)\T) = J Da{x)S[a] W T [a] (6) 

We have introduced here the weight function Wt [a] which should be thought 
of as the square of the target wave function and thus gives the probability 
density to find a particular configuration of a in the wave function of the 
target. There is a subtlety involved here, and the identification of W with 
the square of the wave function is not entirely correct [18]. But it is good 
enough for all we are going to do it these lectures. We can also redefine W 
in such a way that the basic variable becomes S. So 

a(l gluon) = J DS{x)S{x)W T [S] (7) 

with the understanding that W is normalized as the probability density 

j DSW T [S] = 1 (8) 

The functionals Wt in equations © and (JJJ) strictly speaking differ by the 
Jacobean of the transformation but for the simplicity of notation I will 
not indicate this explicitly. 

Consider now an arbitrary projectile whose wave function in the gluon 
Fock space can be written as 

\P) = *[a] a (x)]\0) (9) 

The gluon creation operator a* depends on the transverse coordinate, and 
also on the longitudinal momentum k + . We assume that all the gluons are 
energetic, and therefore all gluon operators which enter eq.@ have longi- 
tudinal momenta above some cutoff A. We will refer to these degrees of 
freedom as "valence". The longitudinal momenta of the gluons are not ex- 
plicitly written down in eq.@ because we assume that as long as those 
gluons are energetic enough they all scatter eikonally, and thus their ampli- 
tude does not depend on the exact value of the momentum. 

Now, since the interaction with the target takes a very short time, all 
gluons scatter independently of each other. Thus the outgoing state simply 
has every gluon creation operator multiplied by the single gluons scattering 
matrix at the corresponding transverse position. 



S\P) = *[S a \x)a] b (x)]\0) 



(10) 
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Given this, the formal expression for the S - matrix of the projectile can 
be written as 

E[S] = (P\S\P) = (0mat(x)MS ab (x)4 b (x)]\0) (11) 
The forward scattering amplitude for the process is thus given by 

s(P on T) = (S)r = J dS E[5] W T [S\ (12) 

Now given these engredients, we ask how does the forward amplitude 
change as we make more energy available for the collision. 

2.2. What of boost? 

We can increase the center of mass energy of the colliding system by 
either boosting the target or the projectile. Of course it does not matter, 
since the S-matrix is Lorentz scalar. The JIMWLK equation, which is 
what we set out to derive, is the equation for the evolution of the target 
probability density function Wt- However it is more convenient to derive it 
by considering boosting the projectile, and then imposing the requirement 
of Lorentz invariance of the amplitude. 

So let's think what happens if we boost the projectile. In what we did so 
far, we assumed that we know the projectile wave function at longitudinal 
momenta greater than some cutoff A. We also tacitly assumed that the 
number of particles below this cutoff is small. Otherwise they would also 
contribute to the S-matrix. So even though we do not know exactly how the 
wave function looks for k + < A, we are not particularly worried. However 
once we start to be interested in boosting, this becomes a very important 
question. The point is that in QCD boosting a state changes the number of 
gluons in it. Let's pause and ask ourselves how does this happen? 

Let's say we boost the state with some rapidity y. We know that one 
effect of this is to scale the longitudinal momenta of all the particles present 
in the wave function by 

k + -► e y k + (13) 

So if we had particles sitting above some cutoff A and the phase space below 
the cutoff was small, then the cutoff moves and the phase space increases. 
The boost drags particles from under the cutoff into the "observable" range 
of momenta. The "unknown" phase space is stretched: 

By itself this is not a big problem. If the number of particles does not 
change, then it does not matter much if their longitudinal momenta have 
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Fig. 5. Growth of the longitudinal phase space under boost. 

increased - in our calculation the scattering amplitude does not depend 
on the value of the longitudinal momentum. The trouble is that in QCD 
things are more complicated. Remember: gluons (or photons) are transverse 
modes of the color field (electromagnetic field). In principle the color field 
also has a longitudinal component - Coulomb field, which is non dynamical. 
However boost does not preserve the longitudinal nature of the field, and 
thus transforms the Coulomb field into live gluons. 

Think as an example of a static charge. It is accompanied by the 
Coulomb electric field 

El{x) = iw (15) 

This is purely longitudinal. Now boost it with a very large rapidity. It's a 
little bit too cumbersome to write down explicit expressions for finite boost, 
so let us formally write it for infinite boost. We know that what you get is 
the electric field which is shrunk to a (5-function in the longitudinal direction, 
while since the total flux is not changed the transverse components decrease 
much slower at large transverse distances. 

47r \r\ s 2irxj_ 

This is the so called Weizsaker- Williams (WW) field. It obviously has trans- 
verse components since d^E 1 ^ 0. So by boosting a benign Coulomb field 
we have created live gluons/photons. We can actually calculate how many 
we have created. We will do it more formally in the next lecture, but for 
now we only have to remember how the gluon creation and annihilation 
operators are related to the electric field 



E\k) = i^4k)[ ai {k) - a\{k)] = iVk+[ ai (k) - a\{k)\ (17) 

where the last line assumes large longitudinal momentum. For the Weizsaker- 
Williams field 

E i (k)=ig^ (18) 



zakopane printed on February 2, 2008 



13 



So we find 

a{k) ~ 9 7k?jk (19) 

and the total number of gluons at a given transverse momentum 

N(k ± ) = J dk+(4(k) ai (k)) = ^J^ = ^y (20) 

That's quite impressive: we have created the number of particles propor- 
tional to the total longitudinal phase space opened up by the boost. 

If we boost by a large y, we create lots of gluons out of the Coulomb 
field. And we do have Coulomb field: we have valence gluons in our wave 
function. They have color charges, and therefore by fiat create Coulomb 
field, and this will turn into new gluons after boost; the new gluons also 
have color charge, and so they carry their own Coulomb field and SO IT 
GOES ON. 

Now, the Coulomb field (or the zero mode, or whatever you want to 
call it) is really everything there is below our cutoff A. To know how many 
gluons will be created by the boost, we therefore need to peak under the 
cutoff. 

The lesson is that in order to find the evolution of S we need to know 
the wave function below A to greater accuracy than simply saying there are 
almost no gluons there. This is our next task. W are now going to address 
the question: given that we know the wave function above A, how do we 
find it below A? 



2.3. The Light Cone Hamiltonian. 

Our aim now is to find the projectile wave function at low longitudinal 
momenta. First of all we have to set up the Hamiltonian calculation. Since 
our projectile moves very fast, the formalism that is most suited for it is the 
Light Cone Quantization. This has been used for many years now in the 
context of the parton model, but also more recently to study the spectrum 
of QCD. A comprehensive review of the formalism is given in |3Uj . 

The essence of this setup is that one quantizes not on the equal time 
surface t = and dynamically evolves states in time, but rather on the 
surface x + = 0. The evolution in x + (which is the analog of the time 
variable) is generated by the operator P~ = H — P%. This is done in the 
light cone gauge A + = 0. Yet another component of the vector potential 
is eliminated, since one of the Maxwell's equations becomes a constraint on 
the light cone: 

d + d + A- + d+diA 1 = (21) 
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which is solved by 

A- = -i^diA* (22) 

This is very convenient, since only the transverse field components are left 
as dynamical degrees of freedom and all non dynamical fields are already 
eliminated. The curios feature of the Light Cone Quantization (LCQ) is that 
it has only half of the number of degrees of freedom compared to the equal 
time quantization. Only positive longitudinal momentum components of 
the field create particles, while negative longitudinal components annihilate 
them. This is understandable, since we want to describe the situation where 
all the particles move right with large energy. A state like this can not 
by definition contain particles with negative momentum, and thus negative 
momentum cannot correspond to particle creation, but only to annihilation. 
The mathematical manifestation of this is that the fields A 1 contain both, 
canonical coordinates and canonical momenta and thus do not commute 
at different values of the longitudinal coordinate x~ . When all is said and 
done, and when the quantization is performed properly, the LCQ formalism 
amounts to the following. 

a) . The "time" coordinate is x + , the "space" coordinates are x~ , X{. 

b) . The basic degrees of freedom are the transverse components of vector 
potential with the commutation relation 

Ai(x-, x), Ai(y~, y)] = % -e{x~ - y~) 5 ab <f< 2 >(x - y) , (23) 

where e(x) is the antisymmetric step function e(x) = ^sign(x). In momen- 
tum space the commutator reads 

Al(k + , x), A{(p+, y)] = ^ 5{k + + p+) S ab & 5^ (x - y) , (24) 

The field A can be represented in the standard way in terms of the creation 
and annihilation operators 

1,00 dk + 1 I A / i+ \ -ik+ 



X Ik x 




'<(*+, x). aj V.y)] = (2ir) <U V> Kk* - P + ) « <2) (* - V 

c) The time evolution of the wave functions in this Hilbert space is gener- 
ated by the following light cone Hamiltonian (LCH)(we switched notations 
from P~ to H for simplicity): 



H = -k[ ^n-(fc+,x)n-(-A;+, X ) + ^(fc + ,x)G^(-fc + ,x) 
2vr Jk+>o L2 4 



(26) 
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where the electric and magnetic pieces are 

II-(x-,x) = ^(£V?+^) a (x-,x) (27) 

x) = ^ a (x-,x)-a^i(^ ) x)- 5 r fcc 4(x-,x)^ c (x-,x) 
d) The physical wave functions are subject to the Gauss' law constraint 

C a (x) = | dx" (D i a + ^) a (x^,x) = 0. (28) 

This constraint is the reflection of the fact that A + = is not a complete 
gauge fixing. It still leaves the residual freedom to perform gauge transfor- 
mations with the gauge function which does not depend on x~ . Hence the 
integral over x~ in the constraint ea.(|28|). 

The standard perturbation theory can be formulated straightforwardly 
in this framework. First of all, one notices that even though the theory is 
interacting, the exact ground state of this Hamiltonian is the Fock space 
vacuum of the operators a(k). This is easy to understand. Since the Hamil- 
tonian commutes with the total momentum operator in the longitudinal 
direction, P + , when written in momentum space each and every term in 
this Hamiltonian has to have vanishing total longitudinal momentum. But 
this immediately means that there can be no terms which contain only cre- 
ation operators of the type (a^) n , since all creation operators have positive 
longitudinal momentum! Thus each term must contain at least one annihi- 
lation operator, and therefore the Fock space vacuum of a is also the vacuum 
of H. 

At any rate the free Hamiltonian is 

H °=LXB?^° ik+Mank+ ^ ) (29) 

and this defines the standard free dispersion relation 

k 2 

2k'k + -k 2 ± = 0, k~ = -rq: (30) 

We are however not interested in the vacuum state. We want to examine 
the state which contains some particles at valence momenta k + > A, and 
in particular to find how this state looks like at momenta below A. There 
is a very natural framework in which to approach this problem, namely 
the Born-Oppenheimer approximation. The free dispersion relation k~ = 
tells us that the frequencies of the high momentum modes are small 
while those of small momentum modes are large. Thus in the framework 
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of the light cone evolution, the high longitudinal momentum modes are 
slow degrees of freedom, and we can study the low momentum mode wave 
function by momentarily freezing the high k + part. This of course is very 
much in line with the picture of the Gluon Cloud which cools and freezes 
under the boost transformation. 

The most important characteristic of the high momentum component of 
the state is its colour charge density 

p a (x ± )= f aJ fe (fc+,x ± )T b X(fc+,x ± ) (31) 

Jk+>A 

since the eikonal 5-matrix depends only on the colour charge density at 
every point in the transverse plane. 

So the strategy is the following. We will forget about the self interaction 
of the high k + modes in the Hamiltonian. We will keep the self interaction 
of low k + modes as is. And additionally we will keep the eikonal interaction 
between the fast and slow modes, namely the interaction between the low 
k + modes and the color charge density generated by the valence modes. 

One can question the applicability of the Born-Oppenheimer approxima- 
tion, since the spectrum of frequencies in the problem is actually continuous, 
so there is no big frequency gap between what we call "fast" and what we 
call "slow". Fortunately however the k + < A bin is dominated by modes 
with momentum much smaller than A. The relevant corrections are propor- 
tional to the total phase space available in the low k + bin, and that goes as 
y = ln(A/Ao). This means that the momenta that dominate the integrals 
are of order k + ~ y/AAo <C A. 

I will not further justify this approximation here, but it turns out to 
be precisely the same eikonal approximation which we used to define the 
scattering matrix. 

All said and done the coupling between the valence and the soft modes 
appears in the Hamiltonian through the modified electric contribution 

n ."<* + ' x ' - =#W' n ^ + - x > + ^wnrf'^ (32) 

+ »ri^/-X +|<A ^ + -P + .x)(-^)^ + ,x)^ 

Now, although p is an operator on the valence Hilbert space, as far as the 
soft fields are concerned, it is a c-number function. Also we will assume, 
consistently with our discussion so far that the number of valence gluons is 
parametrically of order one, and thus also 



p a ~ O(l) 



(33) 
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This makes life much easier, as the interaction term is multiplied by the 
QCD coupling constant. So we can apply perturbation theory. Let us 
therefore calculate the lowest order perturbative correction to the vacuum 
wave function of the soft gluons in the presence of valence charge density. 
As already explained, to zeroth order the vacuum of the LCH is simply the 
Fock space vacuum of a. To first order the correction to the vacuum is given 
by 

(vacuum) = |0) + S^)^M (34 ) 

where the summation is over the eigenstates of the zeroth order Hamiltonian, 
and Ei denotes the zeroth order energy of the respective eigenstate. What 
makes the calculation simple, is that to first order only the interaction term 
that involves the valence charge gives a contribution. The reason is that, as 
we have noted before the Fock space vacuum is actually annihilated by all 
terms of the LCH, as they contain at least one annihilation operator. This 
is not the case however for the soft-valence coupling 



SH 



I tr|^^ [4 a (k + ,k ±)+ a?(k+,k ± )} P a (-k,) (35) 



This Hamiltonian creates only one particle state from the vacuum: 

(k ± ,k + ,a,i\SH\Q) = -J^p°(- k± ) (36) 

With Ei determined by the free dispersion relation we can write the soft 
gluon vacuum state to first order in the coupling as 

Ivacuum) = jl + ij d 2 xb?(x) -^^ Q | a (fc+ *)} |0> (37) 

where the "classical" field b{ is precisely the Weizsaker- Williams field of the 
colour charge density p a 

h ^ x) = hi d2y ^w^ y) - (38) 

Actually one has to be a little bit more careful. We have found the correct 
wave function, but it is not properly normalized. The normalization strictly 
speaking is the second order correction in g, but it has to be restored if we 
want to have unitary evolution. This is easily done, by simply subtracting 
from the coefficient of the Fock vacuum term the total probability to have 
a one gluon state in the wave function. 
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+ i 



I La ^WW~^ {k+ ' x) ) |0) (39) 



In fact, as will become obvious later, this is the only 0{a s ) correction to the 
wave function which gives contribution of 0(a s ) to the evolution equation. 

So we found what we were looking for - the Weizsacker- Williams field 
hiding behind the cutoff. Boosting the state we will pull out this field from 
underneath the cutoff and into the wide open world. Thus starting at the 
initial rapidity Yq with some valence state with momenta above A, call it 
\P)y = \P, k + > A) and boosting it by the boost parameter 5Y we obtain 



1 



P)y +sy - 
I d 2 xb^x)J A 



1 r(l+SY)A fjh.-t r- 

hi /**(*(*)*(*)) 



(40) 



(1+5Y)A dk + 



a\ a (k + ,x) \ \P,k + > (1 + 5Y)A) 



7rV2|fc+|i/2 



We have boosted all the existing gluons to higher momentum, plus have 
produced new gluons all over the newly opened phase space. It is straight- 
forward to see, that the number of newly produced gluons is simply propor- 
tional to the phase space and to the intensity of the Weizsacker- Williams 
field 

dk + (4 a (k + ,x ± )aUk + ,x ± )) = —b»(x)b?(x) (41) 

7T 

2.4. Forty two 

Given that we know now the evolution of the wave function, we can 
find the evolution of the S'-matrix. We go a little back to recall the eikonal 
expression for the S'-matrix 



s(P on T) = (E) T = J dS E[S] W T [S] 

m ^ (P\S\P) = m*K(x)MS ab (x)4\x)]\0) 

The evolved wave function still scatters eikonally. 

r(l+5Y)A dk + f 



(42) 



S\P)y +sy 



1 



-f 



k+ 



I 



d 2 x(b«[Sp,x]b*[S P ,x]) 



(43) 
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+i J d 2 xb?[Sp,x] J^ 1+5Y)A 1 1/2 S ab (x)af(k + , x) | x 

£|p,fc + > (i + <yr)A) 

Here the charge density that creates the Weizsacker- Williams field in the 
outgoing wave function is also rotated by the same one gluon scattering 
matrix 

bt[ X ,S P ] = £jdh^iS^p\y). (44) 

Since the initial projectile wave function does not contain soft gluons, 
the calculation of the averages of the soft gluon operators is straightforward. 
The S-matrix of the boosted projectile is therefore (to first order in SY): 

Zy+syIS] = Ey[S] (45) 

SY f 
' <0|**[<] d 2 z[b?([z,p})b?[z,p] +b?[z,Sp}b?[z,Sp] 



2vr 

-2 bf[z, p] S ba (z) b\[z, Sp}] nSabaf 1 ] |0> 

Our aim is actually to represent this expression in terms of functional 
derivatives of E with respect to S. This involves some algebra, but it is quite 
cute. First we note that the charge density operator, when acting on the 
wave function simply rotates the gluon creation operator at the transverse 
position in question. We can therefore write 

P a (xMs4]\o) = t&?(x)-£— nsaim 

Sal ( x ) 

= -t x {s{x)T*^^}9[Sal]\0) ; (46) 

The second equality is simply the reflection of the fact that in an expression 
of the form Sa, the left multiplication of a is equivalent to the right mul- 
tiplication of S. So the charge density operator acts like the right rotation 
on the matrix S 

p a (x){S bc (x)} = -[S{x)T a ] bc (47) 

Since our wave function is assumed to be dilute (does not contain more than 
one gluon at each point in space) we can immediately write 

p<WW<l>[SaI]|0) = (48) 
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[S(x)p(x)] a ^[Sa\]\0) = -tr\T a S(x)-J—\nSa\}\0) ; 



It is also easy to see that the operator Sp acts as the left rotation: 
p b (y) [S(xMx)] a *[s4}\0) 

= tr tr { s{y)Tb Jw) } *p*> 

etc. (49) 
We can thus immediately write 

^E[S] = X [S,-fg]V[S\ (50) 
where the kernel x is given by 

^JIMWLK _ a s f Cg ~ x )i{ z — y)i / K1 s 
X ~ ^733 I 7Z 1333 3a x ^ oi J 



2 7T 2 Ja;^^ (z - x) 2 (z - 7/) 2 

J a L {x)Jt{y) + J a R (x)J a R (y) ~ 2r L (y)J b R (x)S ba (z) 
with the right- and left rotation operators are defined as 

J?(x) = {S(x)J R (x)]°, (52) 

This is the JIMWLK equation. Actually, those who have seen JIMWLK 
before will probably not quite recognize it in this formula. For starters, 
JIMWLK is usually written as the evolution equation for the target proba- 
bility density function Wt- This is easy to do using Lorentz invariance of the 
forward scattering amplitude. Recall that the forward scattering amplitude 
is given by the equation (|43j) . Let us restore the rapidity labels on various 
quantities. The S matrix is defined at rapidity of interest Y. But we can 
think of this rapidity as distributed between the target and the projectile. 
Let's say that the target has been boosted to rapidity Yq, and the projectile 
to rapidity Y — Yq, so that ea.l|43|) properly reads 

s Y (P on T) = (S y " y ») T(yo) = J dS Z Y - Yo [S] W Yo [S] (53) 

Lorentz invariance requires that s does not depend on Yq. Requiring that 
the derivative of s with respect to Yq vanishes, we find immediately 

A W T [S] = X JIMWLK [S, ±] W T [S] . (54) 
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This is still not the form in which JIMWLK is usually written down. To 
get it into that form, recall that the eikonal expression of the single gluon 
5-matrix is given by 



S(x) = Pexp\i dx-a a (x,x-)T a \ (55) 




with P denoting the path ordering. When acting on any function of S, the 
derivative with respect to a(x, Y) is identical to the action of Jr{x): 



5a e (x,Y) 1 1 5a e (x,Y) 5S ab (x) 



Thus we can write ea . 1)54(1 in the familiar form 

jimwlk _ a s f {z — x)i{z — y)i 



X 27r 2 J x , y , z (z-x) 2 (z-y) 2 { ' 

S\x)S(y) - SHx)S(z) - SHz)S(y)] ab 



5a a {x,Y) L v v w v ; \ j v / Sa b (y,Y) 

Let me finish this lecture with some remarks regarding the JIMWLK 
equation. First, as is clear from the above derivation, the introduction of 
the additional coordinate x~ and the path ordering in eq.© is not really 
necessary. It merely serves to mimic the action of iJ^ on a function of S as 
in ea. ((56() . In the framework of JIMWLK equation it has no independent 
meaning, since none of the physical observables depend on x~ . JIMWLK 
lives happily without the additional coordinate and is defined exclusively 
in terms of the unitary matrix S and its functional derivative 5/5S. I have 
only written it here in this form for those who have seen it before and are 
used to this particular notation. My personal preference is clearly not to 
introduce unnecessary variables, so I prefer the much simpler form ea. 1)52(1. 

The second remark is about the applicability of this equation. We have 
derived it here assuming the projectile is small and deriving the evolution 
of the projectile wave function. On the other hand the original derivation 
was given directly for the target wave function, assuming explicitly that 
the target is dense, that is that the average of the single gluon S matrix 
is close to zero [3]. How come we get the same equations under seemingly 
unrelated conditions? The answer is that these two conditions - diluteness 
of the projectile and dense nature of the target are actually inextricably 
related by the form of the 5-matrix ea. ()43(l . Mathematically diluteness 
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of the projectile is simply the statement that the projectile averaged S- 
matrix S has only low order terms in the Taylor expansion in powers of 
S. On the other hand, if the target is dense it means that the probability 
density functional Wt is strongly peaked around S = 0. If that is the case, 
clearly ea . (|43|) will again be dominated by the first few terms in the Taylor 
expansion of E. Thus the diluteness of projectile, and the dense nature of 
the target in ea, ()43|) go hand in hand. 

Finally the last remark is related to the operatorial nature of S. As 
I have mentioned before, the single gluon scattering matrix is in fact an 
operator in the Hilbert space of the target. Nevertheless we have allowed 
ourselves to treat it as a c- number field neglecting possible commutators. 
The reason we can do this is because the target is dense. What this means, 
is that the charge density of the target is large. Parametrically JIMWLK 
equation corresponds to the target charge densities p ~ In this regime 
the commutator [p a , p b ] = iT^ b p° is suppressed by a s relative to the product 
of p's and can therefore be neglected. This makes a's also classical and 
together with them the matrices S themselves. 

3. JIMWLK and friends 

Let me now discuss some conspicuous limits of the JIMWLK equation 
and some properties of its solutions. 

3.1. The dipole model 

Since its introduction by Al Mueller [27], the dipole degrees of freedom 
have become very popular in the discussion of high energy scattering. The 
formal justification for this model is the large iV c limit. In this limit, in 
any diagram a gluon line can be represented as a double line of a fictitious 
quark and antiquark. It turns out that the diagrams in which these fictitious 
quarks and antiquarks form singlet states - dipoles, are leading in the large 
N c limit. The vertex of emission of the WW gluon then is simply represented 
as the splitting of one dipole into two. 




y 



Fig. 6. In the large N c limit emission of one gluon off a dipole is equivalent to the 
splitting of the dipole into two. 
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The simplicity of the dipole picture, is that there is no color degree of 
freedom anymore, so that the S-matrix of a dipole is a single number rather 
than an SU(N) matrix. Since the color was the only degree of freedom in 
which the eikonal S'-matrix was not diagonal, the S matrix in the dipole 
model is in fact diagonal. We can derive the dipole model limit of the 
JIMWLK equation in a very simple way. The S matrix for a dipole is 
related to the single gluon S'-matrix as: 

s(x,y) = ^Tr(sUx)S F (y)) (58) 

where now Sp is the fundamental representation of the single gluon scat- 
tering matrix (which is the scattering matrix of the fictitious quark which 
constitutes the dipole). We now assume that the projectile and the target 
consist entirely of dipoles, that is that the probability density functional 
W and the projectile scattering matrix £ depends only on this particular 
combination W(s), S(s). It is then straightforward to see that the action 
of the left and right color rotations in the leading N c approximation reduce 
to the following 



J£(x)Jl(y)Z[s] 



J b L {x)r R {y)ns] 



Here r a are the generators of SU{N) group in the fundamental representa- 
tion. We have also dropped the subscript F for brevity, but all the matrices 
in ea.(|59|) are in the fundamental representation. To simplify this further 



± tv[S\u)S(v)r a ] [8(v -y)- 5(u - y)] 

JM C ds{u,v) 



±tr[S\u)T a S(v)] [S(v-y)-S(u-y)] 
Jy c ds[u,v) 

I- tr[S\u)S(v)T a T a ] [5(v -y)- 5{u - y)] 



N, 



x [5(v — x) — 5(u — x)} 



5s(u, v] 



+ 0l K 



l- tv[S\u)S(v)T a T a ] [5(v -y)- 5(u - y)] 



N, 



x [5{v — x) — 5{u — x)] 



5s(u, v) 



° .N? 



I- tr[5 t (u)r ft 5(u)r ] [S(v - y) - 5{u - y)] 



x [S(v — x) — 5(u — x)} 



(59) 
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we need to use the completeness relation 



1 



r a/3 r 7<5 



(60) 



and the identity that relates adjoint and fundamental representations of a 
unitary matrix 



Sf{z) = 2Tr r a S F {z)T b S F {z 



(61) 



Working out the color algebra we find that the evolution operator ^ JIMWLK 
in the large N c limit reduces to 



v JIMWLKro\ 1 



Xdipole [ 



a s 

2^ 



(x - yf 



x,y,z 



(x - zf (z - yf 



x [s(x, y) + s(x,z) s(y,z) 
The large iV c limit evolution 



dY 

can be also written as 
ds(x,y) 



dY 



a, Id z 



Xdi P oie[s] S[s]; 

(x - y) 2 



5s(x,y) 
dW[s] 



dY 



Xdipole[s] W[s] 



(62) 



(63) 



(x — z) 2 (y — z) 2 



[s(x,y) - s(x,z)s(y,z)] (64) 



where the operator s is introduced to indicate that this equation generates 
the hierarchy of equations for the target averages of the powers of s(x,y) - 
multi dipole scattering amplitudes. In particular the single dipole amplitude 
satisfies 



d{s(x,y)) 
dY 



a. I d z- 



(x - yf 



(x — zf (y — zf 
x K s (x,y)} ~ (s(x,z)s{y,z)}} 



(65) 



The forward amplitude is related to the scattering probability by 

s(x, y) = 1 — N(x, y) (66) 

Thus eq. Q65JI is the equation for the evolution of the single dipole scattering 
probability on the target 



d(N(x,y)} 
dY 



a s d z 



(x - y f 



(x — z) 2 (y — zf 



(67) 



x [(N(x, z) + N(y, z) - N(x, y) - N(x, z)N(y, z))] 
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In this way of writing it is obvious that the meaning of the last term is 
simply the probability of simultaneous scattering of both dipoles off the 
target. If we assume that the areas of the target on which the dipoles 
scatter are uncorrelated, we can forget about the target averages and write 
down a closed equation for the single dipole probability 

dNix.y) ( (x - y) 2 , , 

-br 1 = ^\ d \ x - z? ly- z y (68) 

x [N(x, z) + N(y, z) - N(x, y) - N(x, z)N(y, z)\ 

This is the Kovchegov equation [7j. Note that it is derived assuming fac- 
torization of the dipole scattering probabilities. This is best valid when the 
target is dense and when the correlation length of the target fields is small, 
smaller than the distance between the two dipoles. Strictly speaking the 
factorization of the nonlinear term is not valid when the target is far from 
saturation, that is N < 1. However in this regime the nonlinear term is not 
important anyway, and so it came to be that the numerical results obtained 
studying the Kovchegov equation j^J so far are very close to the results 
obtained from the whole JIMWLK hierarchy |32j . 



3.2. The BFKL limit 

Omitting the nonlinear term iV 2 in ()65|) we recover the BFKL equation. 
It is of course always nice to see a familiar face, however on the negative side, 
the BFKL equation is known to be problematic. It violates unitary in more 
than one way. That is to say it violates both the unitarity of the scattering 
probability (the probability can not exceed one) and the unitarity of the 
total cross section (the cross section can not grow faster than the Froissart 
bound a total < -^y 2 in any massive theory). 

The first violation is due to the fact that the solution of the BFKL 
equation for N grows without bound. BFKL is a linear equation and thus 
if the kernel has a positive eigenvalue clearly the probability grows. Indeed 
the BFKL solution at asymptotically large rapidities is 26 

m 2 l n 2 166! 

N(x, b) ~ exp{u;Y - In — (69) 

xqx a z Y 

Here u = Aln2N c a s /n, a 2 = U({3)N c a s /iT and ((n) being the Riemann 
zeta function, with x -the size of the dipole, and b -the impact parameter. 
This unitarity problem is cured by JIMWLK. This is easiest to see on the 
level of Kovchegov equation. Clearly in the presence of the nonlinear term 
the point N(x) = 1 is the fixed point of the equation. If one starts with 
small N, it cannot grow beyond unity. 
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In fact, the evolution works in the way quite consistent with our initial 
expectations about the generation of a typical intrinsic momentum scale. 
In particular, at every rapidity Y, there is a dipole size x s (Y) such that for 
the dipoles of larger sizes the target is black (the scattering probability is 
close to unity), while for smaller size dipoles it is basically transparent. The 
inverse of this scale is the saturation momentum Q S (Y). From the point of 
view of the target evolution, this is the scale on which the target is grainy - 
you probe it at smaller distances, you are liable to see very little since the 
fields are correlated and the phase rotation of q and q are the same, while at 
larger separation the phases rotate randomly and the scattering is almost 
unavoidable. I will discuss properties of Q s some more in a short while. 

However there is also another problem, and that is that the BFKL prob- 
ability amplitude grows large at large impact parameters. In particular in 
the BFKL solution the probability is of order one at impact parameters all 
the way to 

b ~ e KujY (70) 

with k somewhat smaller than one, but not significantly so. As long as N 
is smaller than one, the nonlinear term in the Kovchegov equation is not 
effective, thus it is clear that in this peripheral region the solution of the 
nonlinear equation is similar to that of BFKL. The total cross section is 
given by 

a total = J d 2 bN (pj „ e 2KwY ( 7 -g 

This blatantly violates the Froissart bound. The qualitative picture is that 
of the growth of the dipole due to the long range nature of the WW fields 
|24| . Since perturbative gluon is massless, this growth in the transverse 
space is very fast, which leads to the violation of the Froissart bound. To 
solve this problem one has to understand the real Infrared scale of QCD - 
how is it that gluons do not travel farther than Aqqd from their parents. 
This problem is outside the scope of B-JIMWLK 2 , and so there is no real 
wonder that it stays unresolved. 

3.3. The saturation momentum and the geometric scaling. 

Not all is doom and gloom though. Even though the total cross section 
at very high energies is nonpertubative, it does not mean one cannot extract 
useful physics from JIMWLK. Recall, the main motivation to consider the 
physics of saturation is to study the "minijet physics" around and about 
the saturation scale. This can be done, since the main mechanism of the 
generation of Q s is local and is unaffected by the nonperturbative Coulomb 
tails. The properties of the solutions of the Kovchegov equation have been 
studied very extensively [3^, but the main properties can be understood 
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Fig. 7. Violation of the Froissart bound by the perturbative nonlinear QCD evolu- 
tion. The transverse size of the hadron grows exponentially with rapidity due to 
emission of massless gluons. If initially the projectile hadron missed the target, it 
will hit it after a long evolution due to this growth of transverse size. 



in a very simple way |25j . Basically the solutions to Kovchegov equation is 
the same as that of BFKL for transverse sizes for which N < 1, while for 
sizes where the BFKL solution violates unitarity, the BK solution is simply 
unity. A good approximation to the BFKL solution at large energy is (one 
have to integrate the b dependence in the formula above) 

N(x) = a^A 2 ) 1 ^ exp{uY} (72) 

where A is the scale that characterizes the target at initial rapidity (not 
the longitudinal cutoff of the previous lecture!), and Ao = 0.38 is a pure 
number. The crossover scale Q s is determined by equating this expression 
to a number of order unity. One finds 

x 2 A 2 = exp{--^-Y} (73) 
1 — Ao 

which determines the saturation momentum 

Q 2 S = -I = A 2 exp{-^— Y} (74) 
xj 1 — Ao 

The interesting feature of this solution is that if now we rewrite the ampli- 
tude in terms of Q S (Y) (still in the regime where it is small) we obtain 

N(x) = {x 2 Q 2 s (Y)) 1 -^ (75) 

Thus there is no separate dependence on x and Y, but only on one particular 
combination of the two quantities. This phenomenon came to be known as 
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geometric scaling. Of course, generally one expects such a behavior for 
any saturated system at transverse momenta around Q s . Here however 
the surprising thing is, that due to the fact that the amplitude is a simple 
power, the scaling persists to values of momentum much larger than Q s . A 
simple estimate of the values of x where the scaling breaks down goes as 
follows. There are corrections to the BFKL solutions at small values of x. 
The leading correction is the appearance in the exponent of the same factor 

ln %y This factor becomes comparable with the terms that we kept in 

ea.(|72jl when ln - = ujY . Reexpressing this in terms of the saturation 
momentum, we find that the corrections are small all the way to 

So it turns out that in the BFKL solution geometric scaling persists to 
dipole sizes very much smaller than 1 jQ s . Although this is a very interesting 
property, one has to be careful when trying to assess its phenomenological 
relevance. It does not appear to be a universal prediction of saturation but 
a specific property of the solution of the BFKL equation. Therefore, even 
though the DIS data clearly exhibits such scaling 331 for small values of x, 
its description in terms of the solution of JIMWLK equation is only as good 
as the BFKL fit to DIS data [SU- 



4. What's wrong with JIMWLK? 

What we have accomplished so far is the derivation of the evolution 
equation in the leading order in QCD coupling constant, in the eikonal 
approximation and assuming dilute projectile and dense target. So quite a 
number of restrictions and any one of them when lifted may change results 
dramatically. 

In particular, we know that higher order corrections to BFKL equation 
are very large. And this remains the main threat even if the nonlinearities 
are accounted for. Will the infrared part of emission kill the whole idea 
of perturbatively probing saturation? We don't know, but we also can do 
nothing about it at this point, except try and analyze these corrections |35| . 

What about the eikonal approximation? How long can we rely on it? 
Here we pretty much know that it will break down eventually, but we can 
estimate when it will happen. Assume that we start with a dipole with 
initially large energy Eq. The saturation momentum of the target is what 
governs the recoil effects which are neglected in the eikonal approximation. 
As the saturation momentum becomes of the order of the initial energy, 



Qs(Y) = E 



(77) 
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the eikonal approximation should be abandoned. Even assuming that the 
saturation momentum does grow exponentially (which is certainly not the 
case when next to leading corrections are taken into account P5J this con- 
dition gives 

Y = -ln-^- (78) 
a s A-qcd 

This is pretty far away, so the question is really academic. There is plenty 
of rapidity to go before we reach this limit, and so we can safely use eikonal 
approximation for years to come. 

Now what about the nature of the projectile-target system? This ques- 
tion is the focus of very intense investigation in the last year. It goes under 
various names, like Pomeron loops, fluctuations or wave function saturation 
effects. And this is what I want to talk about in the rest of these lectures. 

So why should we bother? Well, JIMWLK is strictly applicable only 
when the target is dense. Amazingly and surprisingly, its formal low density 
limit reduces to the BFKL equation, which is the correct equation at low 
densities. But what we really want to understand is the transition between 
the dilute and the dense regime. And in this region, for sure things are 
missing from JIMWLK. Let us be faithful to the projectile wave function 
picture, and look at the question from this point of view. Here we have made 
the approximation that the density is small. For example our expression for 
the WW field was first order in the color charge density 

k oc ^p (79) 

This approximation is responsible for the very fast exponential growth of 
the projectile density (and so the saturation momentum). Remember, we 
have found that the number of gluons is proportional to b 2 l We can also 
easily calculate the change in the color charge density itself 

^ = balance + 4^^ ( 80 ) 

In our wave function 

(Sp a (x)Sp a (y)) « f { ( X ~% {y ~% a (u)p a {v) (81) 
Ju,v [x - v) z (y - v) z 

Roughly speaking we have a linear homogeneous equation for the evolution 
of the projectile density (we can be more precise and see that the charge 
density correlator actually satisfies the BFKL equation) 
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The solution to which is clearly an exponential in Y 



p 2 (Y)=ex V {KY}p 2 



(83) 



The exponential growth of the projectile density is the origin of the expo- 
nential growth of the saturation momentum Q s . The physics of this is very 
simple - the number of gluons created by evolution is proportional to the 
number of gluons (color charges) already existing in the state. That's be- 
cause every gluon in the wave function emits independently of the presence 
of other gluons. But of course as a result of such evolution the projectile 
wave function itself becomes dense, at which point independent emission 
approximation is not valid anymore. In fact two things will happen when p 
becomes large, and they combine to give an interesting picture. First, there 
will be nonlinear effects in the emission. Simplistically, we can think of 
this by just assuming that the emissions are happening not in the empty 
space, but in the presence of a non vanishing background field (the WW 
field itself). The emission vertex then should be modified to something like 



where D is the covariant derivative in the 6-background, D = d — b. This 
has the effect that when the field b is large (or better to say for momenta 
smaller than the scale determined by b) the emitted field is not proportional 
to p anymore, but is simply of order unity 



There is also another effect, that I will call "bleaching of color". When the 
projectile is dilute, all gluons are emitted into points in space where there 
were no other gluons previously. This means that the charge density at 
these points monotonically grows. But suppose we are now in the dense 
regime, so that a gluon is emitted with large probability into a point where 
there already exists a large color charge. 

What happens to the color at this point? Well, the color Casimir there 
is as likely to decrease as it is likely to increase. The color charge density 
will random walk! In fact, coupled with our expectation that the emission 
probability of the gluon is constant and independent of density, we expect 
genuine random walk process. 



(84) 




p ~ independent of p 



(85) 



p 2 {Y)~p 2 + KY 



(86) 



Much, much slower than exponential! This clearly must have a large ef- 
fect on physics. We do not understand yet very well what quantities will 
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change significantly and which will not. For example Q s may not change 
too much, since it is mainly sensitive to the range of momenta where p is 
borderline small. But it certainly will effect the structure of final states at 
momenta smaller than Q s and may well also affect the properties of the 
gluon distribution above Q s . 

At any rate looks like fun physics to investigate and understand. 

Incidentally, the bleaching of color, I believe is the same as Pomeron 
loops. The gluon is emitted, but is not actually seen, as its color is bleached 
by the gluon emitted in the next rapidity step. So effectively this is a loop. 

To include the saturation effects in the projectile wave function is there- 
fore the same as including Pomeron loops in the high energy evolution. Do 
we know how to do that? Not so far, although there have been a lot of 
progress in this direction recently. I will tell you a little about the results 
obtained so far, and you will hear more about it from Dionysis' talk later 
onplfij. 



4-1- The simpleminded fix. 

Let's go back to our LCH and see if we can do something when the color 
charge is not too small. 



2vr Jk+>o 
n-(A: + ,x) 

1 



in-(fc + ,x)n-(-A; + ,x) + iG^(fc + ,x)Gjf(-fc+,x) 



-1 



-Me 



■a i nj l (fc + ,x) + 



i 



-i(k+ + ie) 



(87) 



+9- 



-i(k+ + ie) 



f, 



abc 



P+\<A 



dp + 



Al(k + -p + , x )(-ip + )Al(p + , x ) 



GW(s-x) = &Ai(x-,x)-cPAl(x-,x)-gf abc Ai(x-,x)Ai(x-,x 
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As before, the Hamiltonian has a linear term, except that now p is 
not small. We usually solve such problems by shifting the field. It's not 
quite straightforward now since p is itself an operator (albeit on a different 
Hilbert space). However let's forget about it for a moment and treat it like 
a c-number. Then we can get rid of the annoying linear term by the unitary 
transformation 

A -► C ] AC (88) 

with the cloud operator 

C = exp \ij p+>K dy± 2 ^ 1/2 [4 a (P + ,V±) + a.i(P + ,y±j\ (89) 

The "classical" field b now solves the full nonlinear Yang-Mills equation: 

dib l a (x±) = p a (x±) , 

8* K(x ± ) - dj bl(x ± ) - gf abc bl(x ± ) bi(x ± ) = . (90) 

This corresponds to the classical solution for the gauge potential of the 
form 

A?(x) = 9(x-)b?(x ± ) (91) 
of the Yang-Mills equations with the source term 

J +a (x~, x±) = 5(x-)p a (x ± ) (92) 

Note that the field b is a two dimensional pure gauge, so it can be written 
as 

6« = % -tT[ T a U^diU} (93) 

It's easy to see that this works, because with this classical field, both II~ 
and classically vanish, and thus classically this certainly minimizes the 
Hamiltonian. It is also clear that for large p this eliminates the "large" 
contribution to the energy which can be as big as 0(l/a s ). In this approx- 
imation the vacuum wave function for the soft gluons is a coherent state 
given by the action of the cloud operator C on the Fock space vacuum. In 
fact, if we take the extra rapidity interval to be small (small boost), we can 
expand the exponential in C to second order in b. Note that the expansion 
parameter here is not a s anymore, as the field is 0(l/g), but the width 
of the rapidity interval. We can then repeat the derivation of the boosted 
state and find 



\V(Y + 5Y)) 



_ 1.SY J d 2 z(b?(z,\p])b?(z),\p}) 

dk+ 
^\k+\ 1 / 2 



(94) 



+ i I d 2 zb?(z, [p])j^ QK 4 a (k+,z) |*(Y)> 
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The only difference with the expression we derived previously is that the 
classical field now is related to the charge density via a nonlinear equation. 
Again we can continue our previous derivation to find that the kernel of the 
evolution kernel now becomes: 

x jimwlk += _ a- f [q(z,J R mz,J R ) + bt(z,J L )bt(z,J L ) (95) 



2vr 2 

-2b?(z, J L )b?(z, J R )S ba (z)] (96) 
Where as before the right- and left rotation operators are 

JSW = -tr {S(X)T^} ; J£ W = -tr {^M^ 

Jl(x) = \S(x)J R (x)f, (97) 

This time b is a fairly complicated differential operator which satisfies the 
equation 

dib z a (x±) = J R (x±) , 

di b> a (x±) - 8j bi(x x ) - gf abc bl(x ± ) b> c (x ± ) = . (98) 

This expression has been derived by Misha Lublinsky and myself |18j.|19j. 
and also recently rederived by Hatta et.al. in an intriguing paper using the 
formalism of effective action It would be nice of course, if this was the 
whole story, but unfortunately it is not. If one does things more carefully 
it turns out that the shift by the classical field mixes positive and negative 
longitudinal momenta. This means that the vacuum of the transformed 
Hamiltonian is not the Fock space vacuum any more, but rather there is a 
Bogoliubov transformation which transforms one into the other. Thus the 
soft vacuum is not simply C\Q) and more work has to be done to find it. 
These corrections, at least on the face of it are not parametrically suppressed 
relative to the coherent transformation. 

So we do not quite know the evolution kernel yet. There are however 
some interesting hints that we can get from these expressions. The propaga- 
tion of the gluon carries an eikonal factor S. On the other hand the classical 
field of the projectile can be expressed in terms of matrices U, which also 
can be written as eikonal-like factors, but this time in terms of the projectile 
charge density, which as we have seen becomes a derivative with respect to 
the phase of the eikonal factor. This is suggestive of some kind of duality, 
just like duality between x and p for a simple harmonic oscillator. Indeed, 
as it turns out, it is possible to prove that Lorentz invariance plus the con- 
dition that the projectile and target wave functions evolve in the same way 
(projectile-target democracy) requires that the evolution be self dual [2*0] . 
In the last lecture I discuss in more detail the derivation of this self duality 
property. 
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4-2. The self duality of high energy evolution 

We start by considering a general expression for the S-matrix of a pro- 
jectile with the wave function \P) on a target with the wave function |T), 
where the total rapidity of the process is Y. The projectile is assumed to 
be moving to the left with total rapidity Y — Yq (and thus has sizable color 
charge density p~), while the target is moving to the right with total ra- 
pidity Yq (and has large p + ). We assume that the projectile and the target 
contain only partons with large k~ and k + momenta respectively: k~ > A~ 
and k + > A + . The eikonal expression for the S'-matrix reads 

S Y = J Dp +a (x,x-) Wl[p + (x~,x)] ££_y [a], (99) 

where S p is the S'-matrix averaged over the projectile wave function 

S p [a] = ( P \ Ve iJdx- jd*xp- a (.x)<* a {x,x-)\ P y ( 100 ) 

(V denotes path ordering with respect to x~), and W [a] is the weight 
function representing the target, which is related to the target wave function 
in the following way: for an arbitrary operator 0[p + ] 

(T\6[p + (x)]\T) = J Dp +a W T [p + (x-,x)} 0[p + (x,x-)] . (101) 

The field a(x) is the A + component of the vector potential in the light 
cone gauge A~~ = 0. This is the natural gauge from the point of view of 
partonic interpretation of the projectile wave function. In the same gauge 
the target charge density p + is related to a through the solution of the 
classical equations of motion 

a a (x,x~)T a = ±(x-y) {S\y,x-) p+ a (y,x-)T a S(y,x-)}, 

S(x,x~) = 7>exp{i[ X dy-T a a a (x,y-)}. (102) 

where T£ c = if abc is the generator of the SU(N) group in the adjoint rep- 
resentation. The unitary matrix S(x, x~ — > oo) as before has the meaning 
of the scattering matrix of a single gluon from the projectile wave function 
on the target |T). 

These formulae require some explanation. First of all, my notations in 
this lecture are a little different, since I have reversed the direction of motion 
of the target and the projectile. More importantly, I have reintroduced 
here the longitudinal coordinate x~ . Recall, that in the discussion of the 



zakopane printed on February 2, 2008 



35 



JIMWLK equation I have stressed that this coordinate is not necessary. 
However this was in the situation where I have explicitly assumed the target 
to be dense. I am not assuming this in the present discussion. The variable 
x~ plays the role of the ordering variable |18j . The quantities p^ are non 
commuting operators in their respective Hilbert spaces. It turns out however 
that one can trade their non commutativity for an extra coordinate. The 
value of this coordinate simply traces the position of the respective operator 
in the expression whose expectation value one is calculating. The details 
do not matter here (those can be found in HS|), but the upshot is that one 
can treat all the charge density operators like classical quantities, if one 
makes them depend on an extra "ordering" variable and calculates their 
correlation functions via eq. ljlOU) . 

One can also define an analog of W T for the wave function of the pro- 
jectile via 

(P\6[p-(x)]\P) = J Dp- a W p [p-]0[p-(x,x-)]. (103) 

With this definition it is straightforward to see that T, p and W p are related 
through a functional Fourier transform. To represent S as a functional inte- 
gral with weight W p we have to order the factors of the charge density pr in 
the expansion of eq. (|1U0|) , and then endow the charge density p~ (x) with an 
additional coordinate t to turn it into a classical variable. This task is made 
easy by the fact that the ordering of p in ea. (|10U|) follows automatically the 
ordering of the coordinate x~ in the path ordered exponential. Since the 
correlators of p(x,ti) with the weight W p depend only on the ordering of 
the coordinates ti and not their values, we can simply set t = x~ . Once we 
have turned the quantum operators p into the classical variables p(x~), the 
path ordering plays no role anymore, and we thus have 

S P (a) = J Dp a W P [p] e ijdx~ JSxp-(x,x-)a-(x,x-)_ ^ 

We now turn to the discussion of the evolution. The evolution to higher 
energy can be achieved by boosting either the projectile or the target. The 
resulting S-matrix should be the same. This is required by the Lorentz 
invariance of the S'-matrix. Consider first boosting the projectile by a small 
rapidity 5Y. This transformation leads to the change of the projectile S'- 
matrix £ of the form 

= x %,-^-]E p [a) (105) 

The explicit form of the kernel x is known in the large density limit (JIMWLK) 
and in the small density limit (KLWMIJ ^H] - which we have not discussed 
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here), but is not known in general. Substituting eq. (|l()5|) into eq. (|99j) we 

have 

^S Y = J Dp+o(x,x-)Wllp+(x-,x)] {x f k^] 

= J Dp +a {x,x-) \x[a,-^]Wl\p + {x-,x)}} S?_ y Ja]. (106) 

Where the second equality follows by integration by parts. We now impose 
the requirement that the S-matrix does not depend on Yq [22j- Since £ in 
eq. {(nnj) depends on the difference of rapidities, eq. (|l()5j) implies 

= -X*[^]EH. (107) 
Requiring that dS/dY$ = we find that W should satisfy 

^W T = x[a ,±.]W T [p + ] (108) 

Thus we have determined the evolution of the target eq. (|1U8|1 by boosting 
the projectile and requiring Lorentz invariance of the S'-matrix, just as we 
did in the previous lecture. On the other hand the extra energy due to boost 
can be deposited in the target rather than in the projectile. How does W T 
change under boost of the target wave function? To answer this question 
we consider the relation between £ and W together with the evolution of 
X. Referring to eas. (|lU4|) and ()105|) it is obvious that multiplication of T, p 
by a is equivalent to acting on W p by the operator —i5/5p, and acting on 
Ti P by 5 /5a is equivalent to multiplying W p by ip. Additionally, the action 
of ip and —id/6p on W p must be in the reverse order to the action of 5/ 5a 
and a on E p . This means that the evolution of the functional W p is given 
by 

^W p = x [-i±ip]W p [p]. (109) 

Although ea. ()109|) refers to the weight functional representing the projectile 
wave function, the target-projectile democracy requires that the functional 
form of the evolution must be the same for W T . Comparing ea. (|lU8|) and 
ea, (jlU9|) we find that the high energy evolution kernel must, as advertised, 
satisfy the selfduality relation 

Jo) = X ^~ i Jp 1 ^ ' ^ U0 ^ 

To stress it again, this mathematical property of selfduality is the re- 
flection of the requirement that the wave functions of the projectile and the 
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target transform in the same way in their respective frames (and gauges) 
under the Lorentz transformation. 

It is interesting that although we do not know the complete kernel %, 
we can check that the known limits of it are consistent with the duality 
ea. (jllU() . In particular, the duality transformation interchanges the large 
and small density limits, and the kernel x is known in both these limits. In 
the large density limit the pertinent expression is the JIMWLK kernel. It 
can be conveniently written as 

Us f {z-x)i{z-y)j 

XJIMWLK = ^ {z _ x?{z _ y? >< (HI) 

5 5 I - 5 2 5 S ab (z) 5 



8a a (x,0) 5a a (y,0) 5a a (x, 1) 5a a (y, 1) 5a a (x,0) K ' 5a b (y, 1) 



with 



S ab (z) 



-r, i f 1 dx-T c a c (x-,z) 



(112) 



The low density limit of the kernel has recently also been calculated recently 
[TH] , It reads 

_ a s f (z-x)i(z-y)j 
{p a {x, 0)p a (y, 0) + p a (x, l)p a (y, 1) - 2 p a (x, 0)R ab (z)p\y, 1)} . 

with 



R ab {z) 



f 1 r!r~T c S 



ab 



(114) 



Obviously the kernel ea.([58|) transforms into the kernel ea. (|113() by the du- 
ality transformation eq.^TOJ). The duality is also present in the dipole limit 
of the evolution as discussed in j^J and [3*7] . 

Although very interesting, this selfduality property is not powerful 
enough to unambiguously detemine the structure of the kernel. For ex- 
ample the BFKL kernel is self dual, but clearly is not complete. The same 
applies to the approximation to the kernel eq. discussed in 1TBJ - see 
for discussion of this point. We hope however that once the kernel is known 
the self duality property will help us to solve the evolution equation. 



5. So long and thanks for all the fish 

So this is where we stand at the moment. We have understood some 
interesting physics, but there surely is more to come. The immediate future 
will probably bring more efforts to understand the Pomeron loops and the 
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construction of the complete evolution kernel in the eikonal approximation. 
This looks like an achievable goal. But then there is of course the question 
of solving these equations, or at least understanding semiquantitatively the 
changes that the Pomeron loops bring into the game. Amazingly, some 
research into the solutions has already been done, even though the equation 
has not been derived yet , ^3] , jSH] • It will be very interesting to see 
whether these result stay relevant in the framework of the full equation. 

I believe there are also other things that can be analyzed in the eikonal 
approach. In this lectures I have focused on calculations of the total hadronic 
cross section. Come to think of it, this is probably the least interesting of 
all the interesting observables, but it is the simplest one - hence a good 
starting point. But there has been work also on less inclusive quantities, 
like calculations of diffr active cross section |4"0"]. one[3J and two gluon emis- 
sion probabilities |42| . |43 j within eikonal approach. I believe this is the area 
where the use of the wave function techniques should be most advantageous. 
The wave function contains much more information than is needed for the 
total cross section. This additional knowledge is necessary to calculate the 
various gluonic observables, like one and two gluon emission probabilities 
and other less inclusive observables. So I hope the techniques discussed in 
these lectures will find more applications in the not-so-distant future and 
will prove to be worthwhile developing further. 

Finally, I would like to thank Krzysztof Golec-Biernat, Michal Prasza- 
lowicz and Andrzej Bialas for the invitation to give these lectures, and all 
the participants of the school for making it such an enjoyable experience. 
I would also like to thank my collaborators over the years and especially 
Urs Wiedemann and Misha Lublinsky who contributed bulk of the ideas 
presented here while also cracking some good jokes in the process. I hope 
there will be more to come. 

For now though, so long and thanks for all the fish. 
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